Todo:

# load libraries and read data
library(tidyverse)
Registered S3 methods overwritten by 'dbplyr':
  method         from
  print.tbl_lazy     
  print.tbl_sql      
-- Attaching packages ----------------------------------------------------------------------------------------------------------- tidyverse 1.3.0 --
v ggplot2 3.3.2     v purrr   0.3.4
v tibble  3.1.0     v dplyr   1.0.2
v tidyr   1.1.2     v stringr 1.4.0
v readr   1.4.0     v forcats 0.5.0
package 㤼㸱tibble㤼㸲 was built under R version 4.0.4-- Conflicts -------------------------------------------------------------------------------------------------------------- tidyverse_conflicts() --
x dplyr::filter() masks stats::filter()
x dplyr::lag()    masks stats::lag()
library(readxl)
library(plotly)
package 㤼㸱plotly㤼㸲 was built under R version 4.0.5
Attaching package: 㤼㸱plotly㤼㸲

The following object is masked from 㤼㸱package:ggplot2㤼㸲:

    last_plot

The following object is masked from 㤼㸱package:stats㤼㸲:

    filter

The following object is masked from 㤼㸱package:graphics㤼㸲:

    layout
ghg_emissions_clean <- read_csv("data/clean_data/ghg_emissions.csv")

-- Column specification ----------------------------------------------------------------------------------------------------------------------------
cols(
  ccp_mapping = col_character(),
  source_name = col_character(),
  pollutant = col_character(),
  year = col_double(),
  value = col_double(),
  units = col_character()
)

Only emissions - > emissions that are greater than 0

ghg_true_emissions <- ghg_emissions_clean %>% 
  filter(year == max(ghg_emissions_clean$year)) %>% 
  filter(value >= 0) %>% 
  mutate(across(where(is.character), ~str_to_title(.)))
# need to split by pollutant and year
ghg_emissions_clean %>% 
  select(ccp_mapping, source_name) %>% 
  filter(str_detect(source_name, paste("^", ccp_mapping, sep = ""))) %>% 
  unique()
# var_names

# get number of separators in source column
n_breaks <- max(str_count(ghg_emissions_clean$source_name, " - "), na.rm = TRUE)

# number of potential children  = number of breaks + 1
# since: "one" - "two" - "three"
n_children <- n_breaks + 1

# create vector and fill with childnames
child_names <- c()

for (i in seq(1:n_children)) {
  child_name <- paste("child_order_", i, sep = "")
  child_names <- c(child_names, child_name)
}

child_names
[1] "child_order_1" "child_order_2" "child_order_3"
# clean source name column and split into child columns
ghg_wide <- ghg_emissions_clean %>% 
  mutate(source_name = str_to_title(
    str_remove(
      source_name,
      paste("^", ccp_mapping, " - ", sep = "")
      ))) %>% 
  separate(source_name, into = child_names, sep = " - ", fill = "right") %>% 
  rename(child_order_0 = ccp_mapping)
ghg_wide_emissions <- ghg_wide %>%
  filter(!value < 0)

ghg_wide_sinks <- ghg_wide %>% 
  filter(value < 0)

get_child_cols <- function(df, additional_vars, standard_vars = c("value", "units")) {
  temp <- names(df)
  remove <- c(additional_vars, standard_vars)
  
  child_cols <- temp [!temp %in% remove]
}


child_cols <- get_child_cols(ghg_wide_emissions, additional_vars)

child_cols
[1] "child_order_0" "child_order_1" "child_order_2" "child_order_3"
previous_children
NULL
create_child_table(ghg_wide_emissions,
                   current_child = "child_order_0",
                   previous_children = NULL,
                   additional_vars = c("pollutant", "year"))
Error: Problem with `mutate()` input `..1`.
x Input `..1` must be a vector, not a `quosure/formula` object.
i Input `..1` is `current_child`.
Run `rlang::last_error()` to see where the error occurred.
previous_children <- NULL

children <- c(previous_children, "child_order_0")

children
[1] "child_order_0"
ghg_wide_emissions %>% 
  group_by(child_order_0, pollutant, year) %>% 
  summarise(value = sum(value), .groups = "drop_last") %>% 
  # id must be created, parents are null - this is the top level
  mutate(id = child_order_0, parent = "") %>% 
  select(id, label = child_order_0, parent, pollutant, year, value)
NA
ghg_wide_emissions %>%
  group_by(child_order_2, child_order_1, child_order_0, pollutant, year) %>% 
  summarise(value = sum(value), .groups = "drop") %>% 
  mutate(id = paste(child_order_0, child_order_1, child_order_2, sep = " - "),
         parent = paste(child_order_0, child_order_1, sep = " - ")) %>% 
  select(id, label = child_order_2, parent, pollutant, year, value)
ghg_wide_emissions %>%
  group_by(child_order_3, child_order_2, child_order_1, child_order_0, pollutant, year) %>% 
  summarise(value = sum(value), .groups = "drop") %>% 
  mutate(id = paste(child_order_0, child_order_1, child_order_2, child_order_3, sep = " - "),
         parent = paste(child_order_0, child_order_1, child_order_2, sep = " - ")) %>% 
  select(id, label = child_order_3, parent, pollutant, year, value)
create_child_table <- function(df, current_child, previous_children, additional_vars) {
  current_child <- quo(current_child)
  previous_children <- quo(previous_children)
  additional_vars <- quo(additional_vars)
  
  df %>% 
    group_by(!!current_child, rev(!!previous_children), !!additional_vars) %>% 
    summarise(value = sum(value), .groups = 'drop') %>% 
    unite(!!previous_children, col = "parent", sep = " - ") %>% 
    unite(c(!!previous_children, !!current_child), sep = " - ") %>% 
    select(id, label = !!current_child, parent, !!additional_vars, value)
}
names(ghg_wide_emissions) %in% "child_order_[0-9]"
[1] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
  
temp <- names(ghg_wide_emissions)
remove <- c(additional_vars, "value", "units")


temp [!temp %in% remove]
[1] "child_order_0" "child_order_1" "child_order_2" "child_order_3"
create_child_table(ghg_wide_emissions, current_child = "child_order_0",
                   previous_children = "", additional_vars = c("pollutant", "year"))
Error: Invalid index: out of bounds
create_child_table(ghg_wide_emissions, current_child = child_order_0,
                   previous_children = NULL, additional_vars = c(pollutant, year))
Error: Problem with `mutate()` input `..1`.
x Input `..1` must be a vector, not a `quosure/formula` object.
i Input `..1` is `current_child`.
Run `rlang::last_error()` to see where the error occurred.
if (is.null(previous_children)) {
  print("bean")
}
[1] "bean"
vector <- c("child_1", "child_2")

paste(paste(vector, collapse = " - "),"child_3", sep = " - ")
[1] "child_1 - child_2 - child_3"

df %>% group_by(current_child, rev(previous_children), additional_vars) %>% summarise(value = sum(value), .groups = “drop”) %>% mutate(id = paste(previous_children, current_child, sep = " - “), parent = paste(previous_children, sep =” - ")) %>% select(id, label = current_child, parent, additional_vars, value)

vector <- c()

paste(paste(vector, sep = " - "),"child_3", sep = " - ")
[1] " - child_3"
vector <- c()

tibble(
  parent = c("bean", "curd", "whey", "sprout", ""),
  people = c("david", "sasha", "john", "smith", "delilah"),
  shapes = c("triangle", "rectangle", "square", "diamond", "")
) %>% 
  unite(c(1,2,3), col = "id", sep = " - ", remove = FALSE)
tibble(
  parent = c("dad", "bean", "", "fava"),
  child = c("")
) %>% 
  select(child) %>% 
  pull()
[1] "" "" "" ""
# get sector totals for parent df (top level)

ghg_sector_totals <- ghg_true_emissions %>% 
  group_by(ccp_mapping) %>% 
  summarise(sector_total = sum(value), .groups = "drop_last")

parent_df <- ghg_sector_totals %>% 
  mutate(parent = "") %>% 
  select(label = ccp_mapping,
         parent,
         value = sector_total)

# the rest of the data

children_df <- ghg_true_emissions %>% 
  filter(year == max(ghg_true_emissions$year)) %>% 
  group_by(source_name, ccp_mapping) %>% 
  summarise(value = sum(value), .groups = "drop_last") %>% 
  select(label = source_name,
         parent = ccp_mapping,
         value)

# combine into one hierarchical df, create id column

emissions_long <- bind_rows(list(parent_df, children_df)) %>% 
  mutate(id = paste(parent, label, sep = " - "), .before = 'label') %>% 
  mutate(id = str_remove(id, "^ - "))

Other method - seperate column

parents <- emissions_wide %>% 
  filter(is.na(second)) %>% 
  select(id, label, parent, value)
first_born <- emissions_wide %>% 
  filter(!is.na(second)) %>% 
  group_by(second, first) %>% 
  summarise(value = sum(value), .groups = "drop_last") %>% 
  mutate(id = paste(first, second, sep = " - ")) %>% 
  ungroup() %>% 
  select(id,
         label = second,
         parent = first,
         value)
second_born <- emissions_wide %>% 
  filter(!is.na(third)) %>% 
  group_by(third, second, first) %>% 
  summarise(value = sum(value), .groups = "drop_last") %>% 
  mutate(id = paste(first, second, third, sep = " - ")) %>% 
  mutate(parent = paste(first, second, sep = " - ")) %>% 
  ungroup() %>% 
  select(id,
         label = third,
         parent,
         value)
third_born <- emissions_wide %>% 
  filter(!is.na(fourth)) %>% 
  group_by(fourth, third, second, first) %>% 
  summarise(value = sum(value), .groups = "drop_last") %>% 
  mutate(id = paste(first, second, third, fourth, sep = " - ")) %>% 
  mutate(parent = paste(first, second, third, sep = " - ")) %>% 
  ungroup() %>% 
  select(id,
         label = fourth,
         parent,
         value)
n_potential_child_layers <- emissions_long %>% 
  select(id) %>% 
  pull() %>% 
  str_count(" - ") %>% 
  max()
# the number of potential child layers is the maximum included in the dataset
n_potential_child_layers <- emissions_long %>% 
  select(id) %>% 
  pull() %>% 
  str_count(" - ") %>% 
  max()
create_hierarchy_df <- function(df, id_col = "id", id_sep = " - ") {
  n_potential_child_layers <- df %>% 
    select(col) %>% 
    pull() %>% 
    str_count(sep) %>% 
    max()
  
  
}
emissions_long %>% 
  data.frame(stringsAsFactors = FALSE) %>% 
  plot_ly(
    ids = ~id,
    labels = ~label,
    parents = ~parent,
    values = ~value,
    type = "sunburst",
    maxdepth = 2,
    insidetextorientation = 'radial',
    marker=list(colorscale='Portland'),
    text = ~units,
    textinfo='label+percent root+percent entry',
    hoverinfo = paste("%{label}: <br>%{value}",'text')
  )
plot_ly(
  labels = c("Eve", "Seth", "Enos", "Noam", "Awan", "Enoch"),
  parents = c("", "Eve", "Seth", "Seth", "Eve", "Awan"),
  values = c(16, 12, 10, 2, 4, 4),
  type = "sunburst",
  branchvalues = "total"
)
tibble(
  labels = c("Eve", "Seth", "Enos", "Noam", "Awan", "Enoch"),
  parents = c("", "Eve", "Seth", "Seth", "Eve", "Awan"),
  values = c(16, 12, 10, 2, 4, 4)
)
fig <- plot_ly(
  labels = cain_ble$labels,
  parents = cain_ble$parents,
  values = cain_ble$values,
  type = 'sunburst'
)

fig
emissions_long$label[1:10]
 [1] "Agriculture"            "Electricity Generation" "Industry"               "Land use"               "Residential"           
 [6] "Services"               "Transport"              "Waste"                  "Accidental fires"       "Accidental fires"      
d <- data.frame(
    ids = c(
    "North America", "Europe", "Australia", "North America - Football", "Soccer",
    "North America - Rugby", "Europe - Football", "Rugby",
    "Europe - American Football","Australia - Football", "Association",
    "Australian Rules", "Autstralia - American Football", "Australia - Rugby",
    "Rugby League", "Rugby Union"
  ),
  labels = c(
    "North<br>America", "Europe", "Australia", "Football", "Soccer", "Rugby",
    "Football", "Rugby", "American<br>Football", "Football", "Association",
    "Australian<br>Rules", "American<br>Football", "Rugby", "Rugby<br>League",
    "Rugby<br>Union"
  ),
  parents = c(
    "", "", "", "North America", "North America", "North America", "Europe",
    "Europe", "Europe","Australia", "Australia - Football", "Australia - Football",
    "Australia - Football", "Australia - Football", "Australia - Rugby",
    "Australia - Rugby"
  ),
  stringsAsFactors = FALSE
)

fig <- plot_ly(d, ids = ~ids, labels = ~labels, parents = ~parents, type = 'sunburst')

d
class(emissions_long)
[1] "tbl_df"     "tbl"        "data.frame"
class(cain_ble)
[1] "tbl_df"     "tbl"        "data.frame"
class(cain_ble)
diff_order <- ghg_emissions_clean %>% 
  dplyr::group_by(ccp_mapping, year, pollutant) %>% 
  dplyr::summarise(value = sum(value), units = units[1], .groups = "keep") %>%
  filter(pollutant == "CO2") %>% 
  filter(year == min(ghg_emissions_clean$year) |
           year == max(ghg_emissions_clean$year)) %>% 
  ungroup() %>% 
  group_by(ccp_mapping) %>% 
  mutate(diff = lag(value) - value) %>% 
  arrange(diff) %>%
  drop_na() %>% 
  select(ccp_mapping) %>% 
  pull()
ghg_emissions_clean %>% 
  dplyr::group_by(ccp_mapping, year, pollutant) %>% 
  dplyr::summarise(value = sum(value), units = units[1], .groups = "keep") %>%
  filter(pollutant == "CO2") %>%
  mutate(ccp_mapping = factor(ccp_mapping, levels = rev(diff_order))) %>% 
  ggplot() +
  aes(x = year, y = value, fill = ccp_mapping) +
  geom_area() +
  facet_wrap(~pollutant) +
  theme_bw()
ghg_emissions_clean %>% 
  group_by(pollutant) %>% 
  summarise(emissions = sum(value))
plotly::ggplotly(
  ghg_emissions_clean %>% 
    filter(year == 2018) %>%
    filter(pollutant %in% c("CO2", "CH4")) %>%
    ggplot() +
    aes(x = factor(ccp_mapping, levels = rev(levels(factor(ccp_mapping)))),
        y = value, fill = source_name,
        text = paste0('</br> Sector: ', ccp_mapping,
                      '</br> Emissions: ', value,
                      '</br> Source Name: ', source_name)) +
    geom_col(position = "stack") +
    theme_bw() +
    theme(legend.position = "none") +
    labs(x = "Sector",
         y = paste0("Emissions (", ghg_emissions_clean$units[1], ")")) +
    coord_flip(),
    tooltip = 'text'
  )
ghg_emissions_data %>% 
  names()
ghg_emissions_data %>% 
  select()
input <- list()

input$col_choice = "national_communication_categories"
ghg_emissions_clean %>%
  group_by_(input$col_choice, "emission_year") %>% 
  summarise(total_ghg_emissions = sum(emissions)) %>% 
  ggplot() +
  aes(x = EmissionYear, y = total_ghg_emissions, group = `National Communication Categories`, colour = `National Communication Categories`) +
  geom_line() +
  geom_point() +
  scale_x_continuous(breaks = seq(1990,2020,5)) +
  theme(legend.position = 0)
ghg_emissions_data %>% 
  distinct(`National Communication Categories`)
ghg_emissions_data %>% 
  filter(EmissionYear != "BaseYear") %>% 
  mutate(EmissionYear = as.numeric(EmissionYear)) %>% 
  group_by(EmissionYear) %>% 
  summarise(total_ghg_emissions = sum(`Emissions (MtCO2e)`)) %>% 
  ggplot() +
  aes(x = EmissionYear, y = total_ghg_emissions) +
  geom_line() +
  geom_point() +
  scale_x_continuous(breaks = seq(1990,2020,5)) +
  ylim(0, 80) +
  theme(legend.position = 0) +
  theme_bw()
ghg_emissions_data %>% 
  distinct(`CCP mapping`)
ghg_emissions_data %>% 
  filter(EmissionYear != "BaseYear") %>% 
  mutate(EmissionYear = as.numeric(EmissionYear)) %>% 
  group_by(`CCP mapping`, EmissionYear) %>% 
  summarise(total_ghg_emissions = sum(`Emissions (MtCO2e)`)) %>% 
  ggplot() +
  aes(x = EmissionYear, y = total_ghg_emissions, group = `CCP mapping`, colour = `CCP mapping`) +
  geom_line() +
  geom_point() +
  scale_x_continuous(breaks = seq(1990,2020,5))
`summarise()` regrouping output by 'CCP mapping' (override with `.groups` argument)

ghg_emissions_data %>% 
  filter(EmissionYear != "BaseYear") %>% 
  mutate(EmissionYear = as.numeric(EmissionYear)) %>% 
  filter(`National Communication Categories` != `CCP mapping`) %>% 
  select(`National Communication Categories`, `CCP mapping`) %>% 
  unique()
emissions_sankey <- emissions_data %>% 
  select(ccp_mapping, source_name, pollutant, emission_year, emissions, units)

category_id,category_name,subcategory_id,subcategory_name,year,emissions,emission

filtered_df <- ghg_emissions_clean %>% 
  select(ccp_mapping, source_name, pollutant, emission_year, emissions, units) %>% 
  filter(pollutant == "CO2") %>% 
  filter(emission_year == "2005")
total_emissions_for_gas <- filtered_df %>% 
  summarise(sum(emissions)) %>% 
  pull()

total_emissions_by_category <- filtered_df %>% 
  select(-pollutant) %>% 
  group_by(ccp_mapping) %>% 
  summarise(cat_sum = sum(emissions), .groups = 'drop_last')
categories <- filtered_df %>% 
  distinct(ccp_mapping) %>% 
  pull()

n_categories <- length(categories)

sources <- filtered_df %>% 
  distinct(source_name) %>% 
  pull()

n_sources <- length(sources)
n_sources
[1] 159
node_names <- c("Total", categories, sources, "Other")

node_names_df <- data.frame("name" = node_names)

total_sankey_tibble <- total_emissions_by_category %>%
  mutate(total = "Total") %>% 
  mutate(total = match(total, node_names) -1) %>% 
  mutate(ccp_mapping = match(ccp_mapping, node_names) -1) %>% 
  select(source = total,
         target = ccp_mapping,
         value = cat_sum)

total_filtered_emissions <- total_sankey_tibble %>% 
  summarise(sum(value)) %>% 
  pull()

other_emissions <- total_emissions_for_gas - total_filtered_emissions

total_other_sankey_tibble <- tibble(
  "source" = c(0),
  "target" = (match("Other", node_names) -1),
  "value" = c(other_emissions)
)


sub_sankey_tibble <- filtered_df %>% 
  select(- c(units, pollutant, emission_year)) %>% 
  mutate(ccp_mapping = match(ccp_mapping, node_names) -1,
         source_name = match(source_name, node_names) -1)

names(sub_sankey_tibble) = c("source", "target", "value")

sankey_tibble <- total_sankey_tibble %>% 
  bind_rows(sub_sankey_tibble) %>% 
  bind_rows(total_other_sankey_tibble)

links_matrix <- data.frame(as.matrix(sankey_tibble, byrow = TRUE, ncols = 3))

# Add a 'group' column to each connection:
links <- links_matrix %>% 
  mutate(group = case_when(
    source == 0 ~ paste("type_", target, sep = ""),
    source!=0 ~ paste("type_", source, sep = "")
  ))

nodes <- node_names_df
# Add a 'group' column to each node.
# All of them in the same group to make them the same colour
nodes$group <- as.factor(c("my_unique_group"))

emissions <- list()

emissions$nodes <- nodes
emissions$links <- links
node_names <- c("Total", categories, sources, "Other")

node_names_df <- data.frame("name" = node_names)

total_sankey_tibble <- total_emissions_by_category %>%
  mutate(total = "Total") %>% 
  mutate(total = match(total, node_names) -1) %>% 
  mutate(ccp_mapping = match(ccp_mapping, node_names) -1) %>% 
  select(source = total,
         target = ccp_mapping,
         value = cat_sum)

total_filtered_emissions <- total_sankey_tibble %>% 
  summarise(sum(value)) %>% 
  pull()

other_emissions <- total_emissions_for_gas - total_filtered_emissions

total_other_sankey_tibble <- tibble(
  "source" = c(0),
  "target" = (match("Other", node_names) -1),
  "value" = c(other_emissions)
)


sub_sankey_tibble <- filtered_df %>% 
  select(- c(units, pollutant, emission_year)) %>% 
  mutate(ccp_mapping = match(ccp_mapping, node_names) -1,
         source_name = match(source_name, node_names) -1)

names(sub_sankey_tibble) = c("source", "target", "value")

sankey_tibble <- total_sankey_tibble %>% 
  bind_rows(sub_sankey_tibble) %>% 
  bind_rows(total_other_sankey_tibble)

links_matrix <- data.frame(as.matrix(sankey_tibble, byrow = TRUE, ncols = 3))

# Add a 'group' column to each connection:
links <- links_matrix %>% 
  mutate(group = case_when(
    source == 0 ~ paste("type_", target, sep = ""),
    source!=0 ~ paste("type_", source, sep = "")
  ))

nodes <- node_names_df
# Add a 'group' column to each node.
# All of them in the same group to make them the same colour
nodes$group <- as.factor(c("my_unique_group"))

emissions <- list()

emissions$nodes <- nodes
emissions$links <- links
total_filtered_emissions <- total_sankey_tibble %>% 
  summarise(sum(value)) %>% 
  pull()

other_emissions <- total_emissions_for_gas - total_filtered_emissions

total_other_sankey_tibble <- tibble(
  "source" = c(0),
  "target" = (match("Other", node_names) -1),
  "value" = c(other_emissions)
)

sub_sankey_tibble <- filtered_tibble %>% 
  select(-category_id, -subcategory_id, -year) %>% 
  mutate(category_name = match(category_name, node_names) -1,
         subcategory_name = match(subcategory_name, node_names) -1)

names(sub_sankey_tibble) = c("source", "target", "value")

sankey_tibble <- total_sankey_tibble %>% 
  bind_rows(sub_sankey_tibble) %>% 
  bind_rows(total_other_sankey_tibble)

links_matrix <- data.frame(as.matrix(sankey_tibble, byrow = TRUE, ncols = 3))

# Add a 'group' column to each connection:
links <- links_matrix %>% 
  mutate(group = case_when(
    source == 0 ~ paste("type_", target, sep = ""),
    source!=0 ~ paste("type_", source, sep = "")
  ))

nodes <- node_names_df
# Add a 'group' column to each node.
# All of them in the same group to make them the same colour
nodes$group <- as.factor(c("my_unique_group"))

emissions <- list()

emissions$nodes <- nodes
emissions$links <- links
make_sankey_dfs <- function(data, userYear, userGas) {
  n_categories <- data %>% 
    distinct(category_name) %>% 
    nrow()
  
  total_emissions_for_gas <- data %>% 
    filter(emission == userGas()) %>% 
    filter(year == userYear()) %>%
    summarise(sum(emissions)) %>% 
    pull()
  
  filtered_tibble <- data %>%
    filter(emission == userGas()) %>% 
    select(-emission) %>% 
    filter(year == userYear()) %>% 
    filter(emissions > userResolution())
  
  total_emissions_by_cat <- filtered_tibble %>%
    group_by(category_name) %>% 
    summarise(cat_sum = sum(emissions), .groups = 'drop_last')
  
  categories <- filtered_tibble %>%
    distinct(category_name) %>% 
    pull()
  
  subcategories <- filtered_tibble %>%
    distinct(subcategory_name) %>% 
    pull()
  
  node_names <- c("Total", categories, subcategories, "Other")
  
  node_names_df <- data.frame("name" = node_names)
  
  total_sankey_tibble <- total_emissions_by_cat %>%
    mutate(total = "Total") %>% 
    mutate(total = match(total, node_names) -1) %>% 
    mutate(category_name = match(category_name, node_names) -1) %>% 
    select(source = total,
           target = category_name,
           value = cat_sum)
  
  total_filtered_emissions <- total_sankey_tibble %>% 
    summarise(sum(value)) %>% 
    pull()
  
  other_emissions <- total_emissions_for_gas - total_filtered_emissions
  
  total_other_sankey_tibble <- tibble(
    "source" = c(0),
    "target" = (match("Other", node_names) -1),
    "value" = c(other_emissions)
  )
  
  sub_sankey_tibble <- filtered_tibble %>% 
    select(-category_id, -subcategory_id, -year) %>% 
    mutate(category_name = match(category_name, node_names) -1,
           subcategory_name = match(subcategory_name, node_names) -1)
  
  names(sub_sankey_tibble) = c("source", "target", "value")
  
  sankey_tibble <- total_sankey_tibble %>% 
    bind_rows(sub_sankey_tibble) %>% 
    bind_rows(total_other_sankey_tibble)
  
  links_matrix <- data.frame(as.matrix(sankey_tibble, byrow = TRUE, ncols = 3))
  
  # Add a 'group' column to each connection:
  links <- links_matrix %>% 
    mutate(group = case_when(
      source == 0 ~ paste("type_", target, sep = ""),
      source!=0 ~ paste("type_", source, sep = "")
    ))
  
  nodes <- node_names_df
  # Add a 'group' column to each node.
  # All of them in the same group to make them the same colour
  nodes$group <- as.factor(c("my_unique_group"))
  
  emissions <- list()
  
  emissions$nodes <- nodes
  emissions$links <- links
  
  return(emissions)
}
make_sankey_dfs(emissions_sankey, userYear = 2005, userGas = "CH4", userResolution = 50)
LS0tDQp0aXRsZTogIlIgTm90ZWJvb2siDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KIyBUb2RvOg0KDQotIGFkZCBjb2x1bW5zIGZvciB0b3RhbHMgZm9yIGVhY2ggcG9sbHV0YW50L3NlY3Rvcg0KDQpgYGB7cn0NCiMgbG9hZCBsaWJyYXJpZXMgYW5kIHJlYWQgZGF0YQ0KbGlicmFyeSh0aWR5dmVyc2UpDQpsaWJyYXJ5KHJlYWR4bCkNCmxpYnJhcnkocGxvdGx5KQ0KDQpnaGdfZW1pc3Npb25zX2NsZWFuIDwtIHJlYWRfY3N2KCJkYXRhL2NsZWFuX2RhdGEvZ2hnX2VtaXNzaW9ucy5jc3YiKQ0KYGBgDQoNCk9ubHkgZW1pc3Npb25zIC0gPiBlbWlzc2lvbnMgdGhhdCBhcmUgZ3JlYXRlciB0aGFuIDANCg0KYGBge3J9DQpnaGdfdHJ1ZV9lbWlzc2lvbnMgPC0gZ2hnX2VtaXNzaW9uc19jbGVhbiAlPiUgDQogIGZpbHRlcih5ZWFyID09IG1heChnaGdfZW1pc3Npb25zX2NsZWFuJHllYXIpKSAlPiUgDQogIGZpbHRlcih2YWx1ZSA+PSAwKSAlPiUgDQogIG11dGF0ZShhY3Jvc3Mod2hlcmUoaXMuY2hhcmFjdGVyKSwgfnN0cl90b190aXRsZSguKSkpDQpgYGANCg0KYGBge3J9DQojIG5lZWQgdG8gc3BsaXQgYnkgcG9sbHV0YW50IGFuZCB5ZWFyDQpnaGdfZW1pc3Npb25zX2NsZWFuICU+JSANCiAgc2VsZWN0KGNjcF9tYXBwaW5nLCBzb3VyY2VfbmFtZSkgJT4lIA0KICBmaWx0ZXIoc3RyX2RldGVjdChzb3VyY2VfbmFtZSwgcGFzdGUoIl4iLCBjY3BfbWFwcGluZywgc2VwID0gIiIpKSkgJT4lIA0KICB1bmlxdWUoKQ0KYGBgDQoNCmBgYHtyfQ0KIyB2YXJfbmFtZXMNCg0KIyBnZXQgbnVtYmVyIG9mIHNlcGFyYXRvcnMgaW4gc291cmNlIGNvbHVtbg0Kbl9icmVha3MgPC0gbWF4KHN0cl9jb3VudChnaGdfZW1pc3Npb25zX2NsZWFuJHNvdXJjZV9uYW1lLCAiIC0gIiksIG5hLnJtID0gVFJVRSkNCg0KIyBudW1iZXIgb2YgcG90ZW50aWFsIGNoaWxkcmVuICA9IG51bWJlciBvZiBicmVha3MgKyAxDQojIHNpbmNlOiAib25lIiAtICJ0d28iIC0gInRocmVlIg0Kbl9jaGlsZHJlbiA8LSBuX2JyZWFrcyArIDENCg0KIyBjcmVhdGUgdmVjdG9yIGFuZCBmaWxsIHdpdGggY2hpbGRuYW1lcw0KY2hpbGRfbmFtZXMgPC0gYygpDQoNCmZvciAoaSBpbiBzZXEoMTpuX2NoaWxkcmVuKSkgew0KICBjaGlsZF9uYW1lIDwtIHBhc3RlKCJjaGlsZF9vcmRlcl8iLCBpLCBzZXAgPSAiIikNCiAgY2hpbGRfbmFtZXMgPC0gYyhjaGlsZF9uYW1lcywgY2hpbGRfbmFtZSkNCn0NCmBgYA0KDQoNCmBgYHtyfQ0KIyBjbGVhbiBzb3VyY2UgbmFtZSBjb2x1bW4gYW5kIHNwbGl0IGludG8gY2hpbGQgY29sdW1ucw0KZ2hnX3dpZGUgPC0gZ2hnX2VtaXNzaW9uc19jbGVhbiAlPiUgDQogIG11dGF0ZShzb3VyY2VfbmFtZSA9IHN0cl90b190aXRsZSgNCiAgICBzdHJfcmVtb3ZlKA0KICAgICAgc291cmNlX25hbWUsDQogICAgICBwYXN0ZSgiXiIsIGNjcF9tYXBwaW5nLCAiIC0gIiwgc2VwID0gIiIpDQogICAgICApKSkgJT4lIA0KICBzZXBhcmF0ZShzb3VyY2VfbmFtZSwgaW50byA9IGNoaWxkX25hbWVzLCBzZXAgPSAiIC0gIiwgZmlsbCA9ICJyaWdodCIpICU+JSANCiAgcmVuYW1lKGNoaWxkX29yZGVyXzAgPSBjY3BfbWFwcGluZykNCmBgYA0KDQpgYGB7cn0NCmdoZ193aWRlX2VtaXNzaW9ucyA8LSBnaGdfd2lkZSAlPiUNCiAgZmlsdGVyKCF2YWx1ZSA8IDApDQoNCmdoZ193aWRlX3NpbmtzIDwtIGdoZ193aWRlICU+JSANCiAgZmlsdGVyKHZhbHVlIDwgMCkNCmBgYA0KDQpgYGB7cn0NCg0KZ2V0X2NoaWxkX2NvbHMgPC0gZnVuY3Rpb24oZGYsIGFkZGl0aW9uYWxfdmFycywgc3RhbmRhcmRfdmFycyA9IGMoInZhbHVlIiwgInVuaXRzIikpIHsNCiAgdGVtcCA8LSBuYW1lcyhkZikNCiAgcmVtb3ZlIDwtIGMoYWRkaXRpb25hbF92YXJzLCBzdGFuZGFyZF92YXJzKQ0KICANCiAgY2hpbGRfY29scyA8LSB0ZW1wIFshdGVtcCAlaW4lIHJlbW92ZV0NCn0NCg0KDQpjaGlsZF9jb2xzIDwtIGdldF9jaGlsZF9jb2xzKGdoZ193aWRlX2VtaXNzaW9ucywgYWRkaXRpb25hbF92YXJzKQ0KDQpjaGlsZF9jb2xzDQpgYGANCg0KYGBge3J9DQpjaGlsZF9jb2xzIDwtIGMoImNoaWxkX29yZGVyXzAiLCBjaGlsZF9uYW1lcykNCg0KcHJldmlvdXNfY2hpbGRyZW4gPC0gTlVMTA0KDQpoaWVyYXJjaHlfZGYgPC0gdGliYmxlKCkNCg0KZm9yIChjaGlsZCBpbiBjaGlsZF9jb2xzKSB7DQogIHByaW50KGNoaWxkKQ0KICBpbnRlcm1lZGlhdGVfdGliYmxlIDwtIGNyZWF0ZV9jaGlsZF90YWJsZShnaGdfd2lkZV9lbWlzc2lvbnMsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGN1cnJlbnRfY2hpbGQgPSBjaGlsZCwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgcHJldmlvdXNfY2hpbGRyZW4gPSBwcmV2aW91c19jaGlsZHJlbiwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgYWRkaXRpb25hbF92YXJzID0gYygicG9sbHV0YW50IiwgInllYXIiKSkNCiAgcHJldmlvdXNfY2hpbGRyZW4gPC0gYyhwcmV2aW91c19jaGlsZHJlbiwgY2hpbGQpDQogIGhpZXJhcmNoeV9kZiAlPiUgDQogICAgYmluZF9yb3dzKGludGVybWVkaWF0ZV90aWJibGUpDQp9DQoNCmhpZXJhcmNoeV9kZg0KYGBgDQoNCmBgYHtyfQ0KY3JlYXRlX2NoaWxkX3RhYmxlKGdoZ193aWRlX2VtaXNzaW9ucywNCiAgICAgICAgICAgICAgICAgICBjdXJyZW50X2NoaWxkID0gImNoaWxkX29yZGVyXzAiLA0KICAgICAgICAgICAgICAgICAgIHByZXZpb3VzX2NoaWxkcmVuID0gTlVMTCwNCiAgICAgICAgICAgICAgICAgICBhZGRpdGlvbmFsX3ZhcnMgPSBjKCJwb2xsdXRhbnQiLCAieWVhciIpKQ0KYGBgDQoNCg0KYGBge3J9DQpwcmV2aW91c19jaGlsZHJlbiA8LSBOVUxMDQoNCmNoaWxkcmVuIDwtIGMocHJldmlvdXNfY2hpbGRyZW4sICJjaGlsZF9vcmRlcl8wIikNCg0KY2hpbGRyZW4NCmBgYA0KDQoNCg0KYGBge3J9DQpnaGdfd2lkZV9lbWlzc2lvbnMgJT4lIA0KICBncm91cF9ieShjaGlsZF9vcmRlcl8wLCBwb2xsdXRhbnQsIHllYXIpICU+JSANCiAgc3VtbWFyaXNlKHZhbHVlID0gc3VtKHZhbHVlKSwgLmdyb3VwcyA9ICJkcm9wX2xhc3QiKSAlPiUgDQogICMgaWQgbXVzdCBiZSBjcmVhdGVkLCBwYXJlbnRzIGFyZSBudWxsIC0gdGhpcyBpcyB0aGUgdG9wIGxldmVsDQogIG11dGF0ZShpZCA9IGNoaWxkX29yZGVyXzAsIHBhcmVudCA9ICIiKSAlPiUgDQogIHNlbGVjdChpZCwgbGFiZWwgPSBjaGlsZF9vcmRlcl8wLCBwYXJlbnQsIHBvbGx1dGFudCwgeWVhciwgdmFsdWUpDQoNCmBgYA0KDQpgYGB7cn0NCmdoZ193aWRlX2VtaXNzaW9ucyAlPiUNCiAgZ3JvdXBfYnkoY2hpbGRfb3JkZXJfMSwgY2hpbGRfb3JkZXJfMCwgcG9sbHV0YW50LCB5ZWFyKSAlPiUgDQogIHN1bW1hcmlzZSh2YWx1ZSA9IHN1bSh2YWx1ZSksIC5ncm91cHMgPSAiZHJvcCIpICU+JSANCiAgbXV0YXRlKGlkID0gcGFzdGUoY2hpbGRfb3JkZXJfMCwgY2hpbGRfb3JkZXJfMSwgc2VwID0gIiAtICIpLA0KICAgICAgICAgcGFyZW50ID0gY2hpbGRfb3JkZXJfMCkgJT4lIA0KICBzZWxlY3QoaWQsIGxhYmVsID0gY2hpbGRfb3JkZXJfMSwgcGFyZW50LCBwb2xsdXRhbnQsIHllYXIsIHZhbHVlKQ0KYGBgDQoNCmBgYHtyfQ0KZ2hnX3dpZGVfZW1pc3Npb25zICU+JQ0KICBncm91cF9ieShjaGlsZF9vcmRlcl8yLCBjaGlsZF9vcmRlcl8xLCBjaGlsZF9vcmRlcl8wLCBwb2xsdXRhbnQsIHllYXIpICU+JSANCiAgc3VtbWFyaXNlKHZhbHVlID0gc3VtKHZhbHVlKSwgLmdyb3VwcyA9ICJkcm9wIikgJT4lIA0KICBtdXRhdGUoaWQgPSBwYXN0ZShjaGlsZF9vcmRlcl8wLCBjaGlsZF9vcmRlcl8xLCBjaGlsZF9vcmRlcl8yLCBzZXAgPSAiIC0gIiksDQogICAgICAgICBwYXJlbnQgPSBwYXN0ZShjaGlsZF9vcmRlcl8wLCBjaGlsZF9vcmRlcl8xLCBzZXAgPSAiIC0gIikpICU+JSANCiAgc2VsZWN0KGlkLCBsYWJlbCA9IGNoaWxkX29yZGVyXzIsIHBhcmVudCwgcG9sbHV0YW50LCB5ZWFyLCB2YWx1ZSkNCmBgYA0KDQpgYGB7cn0NCmdoZ193aWRlX2VtaXNzaW9ucyAlPiUNCiAgZ3JvdXBfYnkoY2hpbGRfb3JkZXJfMywgY2hpbGRfb3JkZXJfMiwgY2hpbGRfb3JkZXJfMSwgY2hpbGRfb3JkZXJfMCwgcG9sbHV0YW50LCB5ZWFyKSAlPiUgDQogIHN1bW1hcmlzZSh2YWx1ZSA9IHN1bSh2YWx1ZSksIC5ncm91cHMgPSAiZHJvcCIpICU+JSANCiAgbXV0YXRlKGlkID0gcGFzdGUoY2hpbGRfb3JkZXJfMCwgY2hpbGRfb3JkZXJfMSwgY2hpbGRfb3JkZXJfMiwgY2hpbGRfb3JkZXJfMywgc2VwID0gIiAtICIpLA0KICAgICAgICAgcGFyZW50ID0gcGFzdGUoY2hpbGRfb3JkZXJfMCwgY2hpbGRfb3JkZXJfMSwgY2hpbGRfb3JkZXJfMiwgc2VwID0gIiAtICIpKSAlPiUgDQogIHNlbGVjdChpZCwgbGFiZWwgPSBjaGlsZF9vcmRlcl8zLCBwYXJlbnQsIHBvbGx1dGFudCwgeWVhciwgdmFsdWUpDQpgYGANCg0KYGBge3J9DQpjcmVhdGVfY2hpbGRfdGFibGUgPC0gZnVuY3Rpb24oZGYsIGN1cnJlbnRfY2hpbGQsIHByZXZpb3VzX2NoaWxkcmVuLCBhZGRpdGlvbmFsX3ZhcnMpIHsNCiAgaGllcmFyY2h5X2RmIDwtIGdoZ193aWRlX2VtaXNzaW9ucyAlPiUgDQogICAgZ3JvdXBfYnlfKC5kb3RzID0gYyhjdXJyZW50X2NoaWxkLA0KICAgICAgICAgICAgICAgICAgICAgICAgcmV2KHByZXZpb3VzX2NoaWxkcmVuKSwNCiAgICAgICAgICAgICAgICAgICAgICAgIGFkZGl0aW9uYWxfdmFycykpICU+JSANCiAgICBzdW1tYXJpc2UodmFsdWUgPSBzdW0odmFsdWUpLCAuZ3JvdXBzID0gJ2Ryb3AnKQ0KICANCiAgaWYgKGlzLm51bGwocHJldmlvdXNfY2hpbGRyZW4pKSB7DQogICAgaGllcmFyY2h5X2RmICU+JQ0KICAgICAgbXV0YXRlKHBhcmVudCA9ICIiKSAlPiUgDQogICAgICBzZWxlY3QoaWQgPSBjdXJyZW50X2NoaWxkLCBsYWJlbCA9IGN1cnJlbnRfY2hpbGQsDQogICAgICAgICAgICAgcGFyZW50LCBhZGRpdGlvbmFsX3ZhcnMsIHZhbHVlKQ0KICB9IGVsc2Ugew0KICAgIGhpZXJhcmNoeV9kZiAlPiUgDQogICAgICB1bml0ZShwcmV2aW91c19jaGlsZHJlbiwgY29sID0gInBhcmVudCIsIHNlcCA9ICIgLSAiLCByZW1vdmUgPSBGQUxTRSkgJT4lIA0KICAgICAgdW5pdGUoYyhwcmV2aW91c19jaGlsZHJlbiwgY3VycmVudF9jaGlsZCksDQogICAgICAgICAgICBjb2wgPSAiaWQiLCBzZXAgPSAiIC0gIiwgcmVtb3ZlID0gRkFMU0UpICU+JSANCiAgICAgIHNlbGVjdChpZCwgbGFiZWwgPSBjdXJyZW50X2NoaWxkLCBwYXJlbnQsIGFkZGl0aW9uYWxfdmFycywgdmFsdWUpDQogIH0NCn0NCmBgYA0KDQoNCmBgYHtyfQ0KbmFtZXMoZ2hnX3dpZGVfZW1pc3Npb25zKSAlaW4lICJjaGlsZF9vcmRlcl9bMC05XSINCiAgDQp0ZW1wIDwtIG5hbWVzKGdoZ193aWRlX2VtaXNzaW9ucykNCnJlbW92ZSA8LSBjKGFkZGl0aW9uYWxfdmFycywgInZhbHVlIiwgInVuaXRzIikNCg0KDQp0ZW1wIFshdGVtcCAlaW4lIHJlbW92ZV0NCmBgYA0KDQpgYGB7cn0NCmhlbHAobWFncml0dHIpDQpzYGBgDQoNCg0KYGBge3J9DQpjcmVhdGVfY2hpbGRfdGFibGUoZ2hnX3dpZGVfZW1pc3Npb25zLCBjdXJyZW50X2NoaWxkID0gImNoaWxkX29yZGVyXzAiLA0KICAgICAgICAgICAgICAgICAgIHByZXZpb3VzX2NoaWxkcmVuID0gIiIsIGFkZGl0aW9uYWxfdmFycyA9IGMoInBvbGx1dGFudCIsICJ5ZWFyIikpDQpgYGANCg0KDQpgYGB7cn0NCmNyZWF0ZV9jaGlsZF90YWJsZShnaGdfd2lkZV9lbWlzc2lvbnMsIGN1cnJlbnRfY2hpbGQgPSBjaGlsZF9vcmRlcl8wLA0KICAgICAgICAgICAgICAgICAgIHByZXZpb3VzX2NoaWxkcmVuID0gTlVMTCwgYWRkaXRpb25hbF92YXJzID0gYyhwb2xsdXRhbnQsIHllYXIpKQ0KYGBgDQoNCmBgYHtyfQ0KY3VycmVudF9jaGlsZCA8LSAiY2hpbGRfb3JkZXJfMSINCnByZXZpb3VzX2NoaWxkcmVuIDwtICJjaGlsZF9vcmRlcl8wIg0KYWRkaXRpb25hbF92YXJzIDwtIGMoInBvbGx1dGFudCIsICJ5ZWFyIikNCg0KaGllcmFyY2h5X2RmIDwtIGdoZ193aWRlX2VtaXNzaW9ucyAlPiUgDQogIGdyb3VwX2J5XyguZG90cyA9IGMoY3VycmVudF9jaGlsZCwgcmV2KHByZXZpb3VzX2NoaWxkcmVuKSwgYWRkaXRpb25hbF92YXJzKSkgJT4lIA0KICBzdW1tYXJpc2UodmFsdWUgPSBzdW0odmFsdWUpLCAuZ3JvdXBzID0gJ2Ryb3AnKQ0KDQppZiAoaXMubnVsbChwcmV2aW91c19jaGlsZHJlbikpIHsNCiAgaGllcmFyY2h5X2RmICU+JQ0KICAgIG11dGF0ZShwYXJlbnQgPSAiIikgJT4lIA0KICAgIHNlbGVjdChpZCA9IGN1cnJlbnRfY2hpbGQsIGxhYmVsID0gY3VycmVudF9jaGlsZCwgcGFyZW50LCBhZGRpdGlvbmFsX3ZhcnMsIHZhbHVlKQ0KfSBlbHNlIHsNCiAgaGllcmFyY2h5X2RmICU+JSANCiAgICB1bml0ZShwcmV2aW91c19jaGlsZHJlbiwgY29sID0gInBhcmVudCIsIHNlcCA9ICIgLSAiLCByZW1vdmUgPSBGQUxTRSkgJT4lIA0KICAgIHVuaXRlKGMocHJldmlvdXNfY2hpbGRyZW4sIGN1cnJlbnRfY2hpbGQpLCBjb2wgPSAiaWQiLCBzZXAgPSAiIC0gIiwgcmVtb3ZlID0gRkFMU0UpICU+JSANCiAgICBzZWxlY3QoaWQsIGxhYmVsID0gY3VycmVudF9jaGlsZCwgcGFyZW50LCBhZGRpdGlvbmFsX3ZhcnMsIHZhbHVlKQ0KfQ0KICANCmBgYA0KDQpgYGB7cn0NCmlmIChpcy5udWxsKHByZXZpb3VzX2NoaWxkcmVuKSkgew0KICBwcmludCgiYmVhbiIpDQp9DQpgYGANCg0KDQoNCg0KZGYgJT4lDQogIGdyb3VwX2J5KGN1cnJlbnRfY2hpbGQsIHJldihwcmV2aW91c19jaGlsZHJlbiksIGFkZGl0aW9uYWxfdmFycykgJT4lIA0KICBzdW1tYXJpc2UodmFsdWUgPSBzdW0odmFsdWUpLCAuZ3JvdXBzID0gImRyb3AiKSAlPiUgDQogIG11dGF0ZShpZCA9IHBhc3RlKHByZXZpb3VzX2NoaWxkcmVuLCBjdXJyZW50X2NoaWxkLCBzZXAgPSAiIC0gIiksDQogICAgICAgICBwYXJlbnQgPSBwYXN0ZShwcmV2aW91c19jaGlsZHJlbiwgc2VwID0gIiAtICIpKSAlPiUgDQogIHNlbGVjdChpZCwgbGFiZWwgPSBjdXJyZW50X2NoaWxkLCBwYXJlbnQsIGFkZGl0aW9uYWxfdmFycywgdmFsdWUpDQoNCmBgYHtyfQ0KdmVjdG9yIDwtIGMoImNoaWxkXzEiLCAiY2hpbGRfMiIpDQoNCnBhc3RlKHBhc3RlKHZlY3RvciwgY29sbGFwc2UgPSAiIC0gIiksImNoaWxkXzMiLCBzZXAgPSAiIC0gIikNCmBgYA0KDQpgYGB7cn0NCnZlY3RvciA8LSBjKCkNCg0KcGFzdGUocGFzdGUodmVjdG9yLCBzZXAgPSAiIC0gIiksImNoaWxkXzMiLCBzZXAgPSAiIC0gIikNCmBgYA0KDQoNCmBgYHtyfQ0KdmVjdG9yIDwtIGMoKQ0KDQp0aWJibGUoDQogIHBhcmVudCA9IGMoImJlYW4iLCAiY3VyZCIsICJ3aGV5IiwgInNwcm91dCIsICIiKSwNCiAgcGVvcGxlID0gYygiZGF2aWQiLCAic2FzaGEiLCAiam9obiIsICJzbWl0aCIsICJkZWxpbGFoIiksDQogIHNoYXBlcyA9IGMoInRyaWFuZ2xlIiwgInJlY3RhbmdsZSIsICJzcXVhcmUiLCAiZGlhbW9uZCIsICIiKQ0KKSAlPiUgDQogIHVuaXRlKGMoMSwyLDMpLCBjb2wgPSAiaWQiLCBzZXAgPSAiIC0gIiwgcmVtb3ZlID0gRkFMU0UpDQpgYGANCg0KDQpgYGB7cn0NCnRpYmJsZSgNCiAgcGFyZW50ID0gYygiZGFkIiwgImJlYW4iLCAiIiwgImZhdmEiKSwNCiAgY2hpbGQgPSBjKCIiKQ0KKSAlPiUgDQogIHNlbGVjdChjaGlsZCkgJT4lIA0KICBwdWxsKCkNCmBgYA0KDQoNCmBgYHtyfQ0KIyBnZXQgc2VjdG9yIHRvdGFscyBmb3IgcGFyZW50IGRmICh0b3AgbGV2ZWwpDQoNCmdoZ19zZWN0b3JfdG90YWxzIDwtIGdoZ190cnVlX2VtaXNzaW9ucyAlPiUgDQogIGdyb3VwX2J5KGNjcF9tYXBwaW5nKSAlPiUgDQogIHN1bW1hcmlzZShzZWN0b3JfdG90YWwgPSBzdW0odmFsdWUpLCAuZ3JvdXBzID0gImRyb3BfbGFzdCIpDQoNCnBhcmVudF9kZiA8LSBnaGdfc2VjdG9yX3RvdGFscyAlPiUgDQogIG11dGF0ZShwYXJlbnQgPSAiIikgJT4lIA0KICBzZWxlY3QobGFiZWwgPSBjY3BfbWFwcGluZywNCiAgICAgICAgIHBhcmVudCwNCiAgICAgICAgIHZhbHVlID0gc2VjdG9yX3RvdGFsKQ0KDQojIHRoZSByZXN0IG9mIHRoZSBkYXRhDQoNCmNoaWxkcmVuX2RmIDwtIGdoZ190cnVlX2VtaXNzaW9ucyAlPiUgDQogIGZpbHRlcih5ZWFyID09IG1heChnaGdfdHJ1ZV9lbWlzc2lvbnMkeWVhcikpICU+JSANCiAgZ3JvdXBfYnkoc291cmNlX25hbWUsIGNjcF9tYXBwaW5nKSAlPiUgDQogIHN1bW1hcmlzZSh2YWx1ZSA9IHN1bSh2YWx1ZSksIC5ncm91cHMgPSAiZHJvcF9sYXN0IikgJT4lIA0KICBzZWxlY3QobGFiZWwgPSBzb3VyY2VfbmFtZSwNCiAgICAgICAgIHBhcmVudCA9IGNjcF9tYXBwaW5nLA0KICAgICAgICAgdmFsdWUpDQoNCiMgY29tYmluZSBpbnRvIG9uZSBoaWVyYXJjaGljYWwgZGYsIGNyZWF0ZSBpZCBjb2x1bW4NCg0KZW1pc3Npb25zX2xvbmcgPC0gYmluZF9yb3dzKGxpc3QocGFyZW50X2RmLCBjaGlsZHJlbl9kZikpICU+JSANCiAgbXV0YXRlKGlkID0gcGFzdGUocGFyZW50LCBsYWJlbCwgc2VwID0gIiAtICIpLCAuYmVmb3JlID0gJ2xhYmVsJykgJT4lIA0KICBtdXRhdGUoaWQgPSBzdHJfcmVtb3ZlKGlkLCAiXiAtICIpKQ0KYGBgDQoNCiMgT3RoZXIgbWV0aG9kIC0gc2VwZXJhdGUgY29sdW1uDQoNCmBgYHtyfQ0KZW1pc3Npb25zX3dpZGUgPC0gZW1pc3Npb25zX2xvbmcgJT4lIA0KICBzZXBhcmF0ZShpZCwgIiAtICIsIGludG8gPSBjKCJmaXJzdCIsICJzZWNvbmQiLCAidGhpcmQiLCAiZm91cnRoIiksDQogICAgICAgICAgIHJlbW92ZSA9IEZBTFNFLCBleHRyYSA9ICJtZXJnZSIsIGZpbGwgPSAicmlnaHQiKQ0KYGBgDQoNCmBgYHtyfQ0KcGFyZW50cyA8LSBlbWlzc2lvbnNfd2lkZSAlPiUgDQogIGZpbHRlcihpcy5uYShzZWNvbmQpKSAlPiUgDQogIHNlbGVjdChpZCwgbGFiZWwsIHBhcmVudCwgdmFsdWUpDQpgYGANCg0KYGBge3J9DQpmaXJzdF9ib3JuIDwtIGVtaXNzaW9uc193aWRlICU+JSANCiAgZmlsdGVyKCFpcy5uYShzZWNvbmQpKSAlPiUgDQogIGdyb3VwX2J5KHNlY29uZCwgZmlyc3QpICU+JSANCiAgc3VtbWFyaXNlKHZhbHVlID0gc3VtKHZhbHVlKSwgLmdyb3VwcyA9ICJkcm9wX2xhc3QiKSAlPiUgDQogIG11dGF0ZShpZCA9IHBhc3RlKGZpcnN0LCBzZWNvbmQsIHNlcCA9ICIgLSAiKSkgJT4lIA0KICB1bmdyb3VwKCkgJT4lIA0KICBzZWxlY3QoaWQsDQogICAgICAgICBsYWJlbCA9IHNlY29uZCwNCiAgICAgICAgIHBhcmVudCA9IGZpcnN0LA0KICAgICAgICAgdmFsdWUpDQpgYGANCg0KYGBge3J9DQpzZWNvbmRfYm9ybiA8LSBlbWlzc2lvbnNfd2lkZSAlPiUgDQogIGZpbHRlcighaXMubmEodGhpcmQpKSAlPiUgDQogIGdyb3VwX2J5KHRoaXJkLCBzZWNvbmQsIGZpcnN0KSAlPiUgDQogIHN1bW1hcmlzZSh2YWx1ZSA9IHN1bSh2YWx1ZSksIC5ncm91cHMgPSAiZHJvcF9sYXN0IikgJT4lIA0KICBtdXRhdGUoaWQgPSBwYXN0ZShmaXJzdCwgc2Vjb25kLCB0aGlyZCwgc2VwID0gIiAtICIpKSAlPiUgDQogIG11dGF0ZShwYXJlbnQgPSBwYXN0ZShmaXJzdCwgc2Vjb25kLCBzZXAgPSAiIC0gIikpICU+JSANCiAgdW5ncm91cCgpICU+JSANCiAgc2VsZWN0KGlkLA0KICAgICAgICAgbGFiZWwgPSB0aGlyZCwNCiAgICAgICAgIHBhcmVudCwNCiAgICAgICAgIHZhbHVlKQ0KYGBgDQoNCmBgYHtyfQ0KdGhpcmRfYm9ybiA8LSBlbWlzc2lvbnNfd2lkZSAlPiUgDQogIGZpbHRlcighaXMubmEoZm91cnRoKSkgJT4lIA0KICBncm91cF9ieShmb3VydGgsIHRoaXJkLCBzZWNvbmQsIGZpcnN0KSAlPiUgDQogIHN1bW1hcmlzZSh2YWx1ZSA9IHN1bSh2YWx1ZSksIC5ncm91cHMgPSAiZHJvcF9sYXN0IikgJT4lIA0KICBtdXRhdGUoaWQgPSBwYXN0ZShmaXJzdCwgc2Vjb25kLCB0aGlyZCwgZm91cnRoLCBzZXAgPSAiIC0gIikpICU+JSANCiAgbXV0YXRlKHBhcmVudCA9IHBhc3RlKGZpcnN0LCBzZWNvbmQsIHRoaXJkLCBzZXAgPSAiIC0gIikpICU+JSANCiAgdW5ncm91cCgpICU+JSANCiAgc2VsZWN0KGlkLA0KICAgICAgICAgbGFiZWwgPSBmb3VydGgsDQogICAgICAgICBwYXJlbnQsDQogICAgICAgICB2YWx1ZSkNCmBgYA0KDQpgYGB7cn0NCiMgdGhlIG51bWJlciBvZiBwb3RlbnRpYWwgY2hpbGQgbGF5ZXJzIGlzIHRoZSBtYXhpbXVtIGluY2x1ZGVkIGluIHRoZSBkYXRhc2V0DQpuX3BvdGVudGlhbF9jaGlsZF9sYXllcnMgPC0gZW1pc3Npb25zX2xvbmcgJT4lIA0KICBzZWxlY3QoaWQpICU+JSANCiAgcHVsbCgpICU+JSANCiAgc3RyX2NvdW50KCIgLSAiKSAlPiUgDQogIG1heCgpDQpgYGANCg0KDQoNCmBgYHtyfQ0KY3JlYXRlX2hpZXJhcmNoeV9kZiA8LSBmdW5jdGlvbihkZiwgaWRfY29sID0gImlkIiwgaWRfc2VwID0gIiAtICIpIHsNCiAgbl9wb3RlbnRpYWxfY2hpbGRfbGF5ZXJzIDwtIGRmICU+JSANCiAgICBzZWxlY3QoY29sKSAlPiUgDQogICAgcHVsbCgpICU+JSANCiAgICBzdHJfY291bnQoc2VwKSAlPiUgDQogICAgbWF4KCkNCiAgDQogIA0KfQ0KYGBgDQoNCg0KDQpgYGB7cn0NCmBgYA0KDQoNCmBgYHtyfQ0KZW1pc3Npb25zX2xvbmcgPC0gYmluZF9yb3dzKGxpc3QocGFyZW50cywgZmlyc3RfYm9ybiwgc2Vjb25kX2Jvcm4sIHRoaXJkX2Jvcm4pKSAlPiUgDQogIG11dGF0ZSh1bml0cyA9IGdoZ19lbWlzc2lvbnNfY2xlYW4kdW5pdHNbMV0pICU+JSANCiAgd3JpdGVfY3N2KCJkYXRhL2NsZWFuX2RhdGEvaGllcmFyY2hpY2FsX2RhdGEuY3N2IikNCmBgYA0KDQoNCg0KYGBge3J9DQplbWlzc2lvbnNfbG9uZyAlPiUNCiAgZGF0YS5mcmFtZShzdHJpbmdzQXNGYWN0b3JzID0gRkFMU0UpICU+JSANCiAgcGxvdF9seSgNCiAgICBpZHMgPSB+aWQsDQogICAgbGFiZWxzID0gfmxhYmVsLA0KICAgIHBhcmVudHMgPSB+cGFyZW50LA0KICAgIHZhbHVlcyA9IH52YWx1ZSwNCiAgICB0eXBlID0gInRyZWVtYXAiLA0KICAgIGJyYW5jaHZhbHVlcyA9ICJ0b3RhbCIsDQogICAgbWF4ZGVwdGggPSAyLA0KICAgIHRleHRpbmZvPSdsYWJlbCtwZXJjZW50IHJvb3QrZW50cnknDQogICkNCmBgYA0KDQoNCmBgYHtyfQ0KZW1pc3Npb25zX2xvbmcgJT4lIA0KICBkYXRhLmZyYW1lKHN0cmluZ3NBc0ZhY3RvcnMgPSBGQUxTRSkgJT4lIA0KICBwbG90X2x5KA0KICAgIGlkcyA9IH5pZCwNCiAgICBsYWJlbHMgPSB+bGFiZWwsDQogICAgcGFyZW50cyA9IH5wYXJlbnQsDQogICAgdmFsdWVzID0gfnZhbHVlLA0KICAgIHR5cGUgPSAic3VuYnVyc3QiLA0KICAgIG1heGRlcHRoID0gMiwNCiAgICBpbnNpZGV0ZXh0b3JpZW50YXRpb24gPSAncmFkaWFsJywNCiAgICBtYXJrZXI9bGlzdChjb2xvcnNjYWxlPSdQb3J0bGFuZCcpLA0KICAgIHRleHQgPSB+dW5pdHMsDQogICAgdGV4dGluZm89J2xhYmVsK3BlcmNlbnQgcm9vdCtwZXJjZW50IGVudHJ5JywNCiAgICBob3ZlcmluZm8gPSBwYXN0ZSgiJXtsYWJlbH06IDxicj4le3ZhbHVlfSIsJ3RleHQnKQ0KICApDQpgYGANCg0KDQoNCmBgYHtyfQ0KcGxvdF9seSgNCiAgbGFiZWxzID0gYygiRXZlIiwgIlNldGgiLCAiRW5vcyIsICJOb2FtIiwgIkF3YW4iLCAiRW5vY2giKSwNCiAgcGFyZW50cyA9IGMoIiIsICJFdmUiLCAiU2V0aCIsICJTZXRoIiwgIkV2ZSIsICJBd2FuIiksDQogIHZhbHVlcyA9IGMoMTYsIDEyLCAxMCwgMiwgNCwgNCksDQogIHR5cGUgPSAic3VuYnVyc3QiLA0KICBicmFuY2h2YWx1ZXMgPSAidG90YWwiDQopDQpgYGANCg0KYGBge3J9DQp0aWJibGUoDQogIGxhYmVscyA9IGMoIkV2ZSIsICJTZXRoIiwgIkVub3MiLCAiTm9hbSIsICJBd2FuIiwgIkVub2NoIiksDQogIHBhcmVudHMgPSBjKCIiLCAiRXZlIiwgIlNldGgiLCAiU2V0aCIsICJFdmUiLCAiQXdhbiIpLA0KICB2YWx1ZXMgPSBjKDE2LCAxMiwgMTAsIDIsIDQsIDQpDQopDQpgYGANCg0KYGBge3J9DQplbWlzc2lvbnNfbG9uZw0KYGBgDQoNCg0KYGBge3J9DQpmaWcgPC0gcGxvdF9seSgNCiAgbGFiZWxzID0gY2Fpbl9ibGUkbGFiZWxzLA0KICBwYXJlbnRzID0gY2Fpbl9ibGUkcGFyZW50cywNCiAgdmFsdWVzID0gY2Fpbl9ibGUkdmFsdWVzLA0KICB0eXBlID0gJ3N1bmJ1cnN0Jw0KKQ0KDQpmaWcNCmBgYA0KDQoNCmBgYHtyfQ0KZW1pc3Npb25zX2xvbmckbGFiZWxbMToxMF0NCmBgYA0KDQoNCmBgYHtyfQ0KZCA8LSBkYXRhLmZyYW1lKA0KICAgIGlkcyA9IGMoDQogICAgIk5vcnRoIEFtZXJpY2EiLCAiRXVyb3BlIiwgIkF1c3RyYWxpYSIsICJOb3J0aCBBbWVyaWNhIC0gRm9vdGJhbGwiLCAiU29jY2VyIiwNCiAgICAiTm9ydGggQW1lcmljYSAtIFJ1Z2J5IiwgIkV1cm9wZSAtIEZvb3RiYWxsIiwgIlJ1Z2J5IiwNCiAgICAiRXVyb3BlIC0gQW1lcmljYW4gRm9vdGJhbGwiLCJBdXN0cmFsaWEgLSBGb290YmFsbCIsICJBc3NvY2lhdGlvbiIsDQogICAgIkF1c3RyYWxpYW4gUnVsZXMiLCAiQXV0c3RyYWxpYSAtIEFtZXJpY2FuIEZvb3RiYWxsIiwgIkF1c3RyYWxpYSAtIFJ1Z2J5IiwNCiAgICAiUnVnYnkgTGVhZ3VlIiwgIlJ1Z2J5IFVuaW9uIg0KICApLA0KICBsYWJlbHMgPSBjKA0KICAgICJOb3J0aDxicj5BbWVyaWNhIiwgIkV1cm9wZSIsICJBdXN0cmFsaWEiLCAiRm9vdGJhbGwiLCAiU29jY2VyIiwgIlJ1Z2J5IiwNCiAgICAiRm9vdGJhbGwiLCAiUnVnYnkiLCAiQW1lcmljYW48YnI+Rm9vdGJhbGwiLCAiRm9vdGJhbGwiLCAiQXNzb2NpYXRpb24iLA0KICAgICJBdXN0cmFsaWFuPGJyPlJ1bGVzIiwgIkFtZXJpY2FuPGJyPkZvb3RiYWxsIiwgIlJ1Z2J5IiwgIlJ1Z2J5PGJyPkxlYWd1ZSIsDQogICAgIlJ1Z2J5PGJyPlVuaW9uIg0KICApLA0KICBwYXJlbnRzID0gYygNCiAgICAiIiwgIiIsICIiLCAiTm9ydGggQW1lcmljYSIsICJOb3J0aCBBbWVyaWNhIiwgIk5vcnRoIEFtZXJpY2EiLCAiRXVyb3BlIiwNCiAgICAiRXVyb3BlIiwgIkV1cm9wZSIsIkF1c3RyYWxpYSIsICJBdXN0cmFsaWEgLSBGb290YmFsbCIsICJBdXN0cmFsaWEgLSBGb290YmFsbCIsDQogICAgIkF1c3RyYWxpYSAtIEZvb3RiYWxsIiwgIkF1c3RyYWxpYSAtIEZvb3RiYWxsIiwgIkF1c3RyYWxpYSAtIFJ1Z2J5IiwNCiAgICAiQXVzdHJhbGlhIC0gUnVnYnkiDQogICksDQogIHN0cmluZ3NBc0ZhY3RvcnMgPSBGQUxTRQ0KKQ0KDQpmaWcgPC0gcGxvdF9seShkLCBpZHMgPSB+aWRzLCBsYWJlbHMgPSB+bGFiZWxzLCBwYXJlbnRzID0gfnBhcmVudHMsIHR5cGUgPSAnc3VuYnVyc3QnKQ0KDQpkDQpgYGANCg0KYGBge3J9DQpjbGFzcyhlbWlzc2lvbnNfbG9uZykNCmBgYA0KYGBge3J9DQpjbGFzcyhjYWluX2JsZSkNCmBgYA0KDQoNCmBgYHtyfQ0KZGlmZl9vcmRlciA8LSBnaGdfZW1pc3Npb25zX2NsZWFuICU+JSANCiAgZHBseXI6Omdyb3VwX2J5KGNjcF9tYXBwaW5nLCB5ZWFyLCBwb2xsdXRhbnQpICU+JSANCiAgZHBseXI6OnN1bW1hcmlzZSh2YWx1ZSA9IHN1bSh2YWx1ZSksIHVuaXRzID0gdW5pdHNbMV0sIC5ncm91cHMgPSAia2VlcCIpICU+JQ0KICBmaWx0ZXIocG9sbHV0YW50ID09ICJDTzIiKSAlPiUgDQogIGZpbHRlcih5ZWFyID09IG1pbihnaGdfZW1pc3Npb25zX2NsZWFuJHllYXIpIHwNCiAgICAgICAgICAgeWVhciA9PSBtYXgoZ2hnX2VtaXNzaW9uc19jbGVhbiR5ZWFyKSkgJT4lIA0KICB1bmdyb3VwKCkgJT4lIA0KICBncm91cF9ieShjY3BfbWFwcGluZykgJT4lIA0KICBtdXRhdGUoZGlmZiA9IGxhZyh2YWx1ZSkgLSB2YWx1ZSkgJT4lIA0KICBhcnJhbmdlKGRpZmYpICU+JQ0KICBkcm9wX25hKCkgJT4lIA0KICBzZWxlY3QoY2NwX21hcHBpbmcpICU+JSANCiAgcHVsbCgpDQpgYGANCg0KDQpgYGB7cn0NCmdoZ19lbWlzc2lvbnNfY2xlYW4gJT4lIA0KICBkcGx5cjo6Z3JvdXBfYnkoY2NwX21hcHBpbmcsIHllYXIsIHBvbGx1dGFudCkgJT4lIA0KICBkcGx5cjo6c3VtbWFyaXNlKHZhbHVlID0gc3VtKHZhbHVlKSwgdW5pdHMgPSB1bml0c1sxXSwgLmdyb3VwcyA9ICJrZWVwIikgJT4lDQogIGZpbHRlcihwb2xsdXRhbnQgPT0gIkNPMiIpICU+JQ0KICBtdXRhdGUoY2NwX21hcHBpbmcgPSBmYWN0b3IoY2NwX21hcHBpbmcsIGxldmVscyA9IHJldihkaWZmX29yZGVyKSkpICU+JSANCiAgZ2dwbG90KCkgKw0KICBhZXMoeCA9IHllYXIsIHkgPSB2YWx1ZSwgZmlsbCA9IGNjcF9tYXBwaW5nKSArDQogIGdlb21fYXJlYSgpICsNCiAgZmFjZXRfd3JhcCh+cG9sbHV0YW50KSArDQogIHRoZW1lX2J3KCkNCmBgYA0KYGBge3J9DQpnaGdfZW1pc3Npb25zX2NsZWFuICU+JSANCiAgZ3JvdXBfYnkocG9sbHV0YW50KSAlPiUgDQogIHN1bW1hcmlzZShlbWlzc2lvbnMgPSBzdW0odmFsdWUpKQ0KYGBgDQoNCg0KYGBge3J9DQpwbG90bHk6OmdncGxvdGx5KA0KICBnaGdfZW1pc3Npb25zX2NsZWFuICU+JSANCiAgICBmaWx0ZXIoeWVhciA9PSAyMDE4KSAlPiUNCiAgICBmaWx0ZXIocG9sbHV0YW50ICVpbiUgYygiQ08yIiwgIkNINCIpKSAlPiUNCiAgICBnZ3Bsb3QoKSArDQogICAgYWVzKHggPSBmYWN0b3IoY2NwX21hcHBpbmcsIGxldmVscyA9IHJldihsZXZlbHMoZmFjdG9yKGNjcF9tYXBwaW5nKSkpKSwNCiAgICAgICAgeSA9IHZhbHVlLCBmaWxsID0gc291cmNlX25hbWUsDQogICAgICAgIHRleHQgPSBwYXN0ZTAoJzwvYnI+IFNlY3RvcjogJywgY2NwX21hcHBpbmcsDQogICAgICAgICAgICAgICAgICAgICAgJzwvYnI+IEVtaXNzaW9uczogJywgdmFsdWUsDQogICAgICAgICAgICAgICAgICAgICAgJzwvYnI+IFNvdXJjZSBOYW1lOiAnLCBzb3VyY2VfbmFtZSkpICsNCiAgICBnZW9tX2NvbChwb3NpdGlvbiA9ICJzdGFjayIpICsNCiAgICB0aGVtZV9idygpICsNCiAgICB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAibm9uZSIpICsNCiAgICBsYWJzKHggPSAiU2VjdG9yIiwNCiAgICAgICAgIHkgPSBwYXN0ZTAoIkVtaXNzaW9ucyAoIiwgZ2hnX2VtaXNzaW9uc19jbGVhbiR1bml0c1sxXSwgIikiKSkgKw0KICAgIGNvb3JkX2ZsaXAoKSwNCiAgICB0b29sdGlwID0gJ3RleHQnDQogICkNCg0KYGBgDQoNCg0KYGBge3J9DQpnaGdfZW1pc3Npb25zX2RhdGEgJT4lIA0KICBuYW1lcygpDQpgYGANCg0KYGBge3J9DQpnaGdfZW1pc3Npb25zX2RhdGEgJT4lIA0KICBzZWxlY3QoKQ0KYGBgDQoNCg0KDQpgYGB7cn0NCmlucHV0IDwtIGxpc3QoKQ0KDQppbnB1dCRjb2xfY2hvaWNlID0gIm5hdGlvbmFsX2NvbW11bmljYXRpb25fY2F0ZWdvcmllcyINCmBgYA0KDQoNCg0KYGBge3J9DQpnaGdfZW1pc3Npb25zX2NsZWFuICU+JQ0KICBncm91cF9ieV8oaW5wdXQkY29sX2Nob2ljZSwgImVtaXNzaW9uX3llYXIiKSAlPiUgDQogIHN1bW1hcmlzZSh0b3RhbF9naGdfZW1pc3Npb25zID0gc3VtKGVtaXNzaW9ucykpICU+JSANCiAgZ2dwbG90KCkgKw0KICBhZXMoeCA9IEVtaXNzaW9uWWVhciwgeSA9IHRvdGFsX2doZ19lbWlzc2lvbnMsIGdyb3VwID0gYE5hdGlvbmFsIENvbW11bmljYXRpb24gQ2F0ZWdvcmllc2AsIGNvbG91ciA9IGBOYXRpb25hbCBDb21tdW5pY2F0aW9uIENhdGVnb3JpZXNgKSArDQogIGdlb21fbGluZSgpICsNCiAgZ2VvbV9wb2ludCgpICsNCiAgc2NhbGVfeF9jb250aW51b3VzKGJyZWFrcyA9IHNlcSgxOTkwLDIwMjAsNSkpICsNCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gMCkNCmBgYA0KDQpgYGB7cn0NCmdoZ19lbWlzc2lvbnNfZGF0YSAlPiUgDQogIGRpc3RpbmN0KGBOYXRpb25hbCBDb21tdW5pY2F0aW9uIENhdGVnb3JpZXNgKQ0KYGBgDQpgYGB7cn0NCmdoZ19lbWlzc2lvbnNfZGF0YSAlPiUgDQogIGZpbHRlcihFbWlzc2lvblllYXIgIT0gIkJhc2VZZWFyIikgJT4lIA0KICBtdXRhdGUoRW1pc3Npb25ZZWFyID0gYXMubnVtZXJpYyhFbWlzc2lvblllYXIpKSAlPiUgDQogIGdyb3VwX2J5KEVtaXNzaW9uWWVhcikgJT4lIA0KICBzdW1tYXJpc2UodG90YWxfZ2hnX2VtaXNzaW9ucyA9IHN1bShgRW1pc3Npb25zIChNdENPMmUpYCkpICU+JSANCiAgZ2dwbG90KCkgKw0KICBhZXMoeCA9IEVtaXNzaW9uWWVhciwgeSA9IHRvdGFsX2doZ19lbWlzc2lvbnMpICsNCiAgZ2VvbV9saW5lKCkgKw0KICBnZW9tX3BvaW50KCkgKw0KICBzY2FsZV94X2NvbnRpbnVvdXMoYnJlYWtzID0gc2VxKDE5OTAsMjAyMCw1KSkgKw0KICB5bGltKDAsIDgwKSArDQogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9IDApICsNCiAgdGhlbWVfYncoKQ0KYGBgDQoNCg0KYGBge3J9DQpnaGdfZW1pc3Npb25zX2RhdGEgJT4lIA0KICBkaXN0aW5jdChgQ0NQIG1hcHBpbmdgKQ0KYGBgDQoNCmBgYHtyfQ0KDQpgYGANCg0KDQpgYGB7cn0NCmdoZ19lbWlzc2lvbnNfZGF0YSAlPiUgDQogIGRpc3RpbmN0KGBOYXRpb25hbCBDb21tdW5pY2F0aW9uIENhdGVnb3JpZXNgKQ0KYGBgDQpgYGB7cn0NCmdoZ19lbWlzc2lvbnNfZGF0YQ0KYGBgDQpgYGB7cn0NCmdoZ19lbWlzc2lvbnNfZGF0YSAlPiUgDQogIGZpbHRlcihFbWlzc2lvblllYXIgIT0gIkJhc2VZZWFyIikgJT4lIA0KICBtdXRhdGUoRW1pc3Npb25ZZWFyID0gYXMubnVtZXJpYyhFbWlzc2lvblllYXIpKSAlPiUgDQogIGdyb3VwX2J5KGBDQ1AgbWFwcGluZ2AsIEVtaXNzaW9uWWVhcikgJT4lIA0KICBzdW1tYXJpc2UodG90YWxfZ2hnX2VtaXNzaW9ucyA9IHN1bShgRW1pc3Npb25zIChNdENPMmUpYCkpICU+JSANCiAgZ2dwbG90KCkgKw0KICBhZXMoeCA9IEVtaXNzaW9uWWVhciwgeSA9IHRvdGFsX2doZ19lbWlzc2lvbnMsIGdyb3VwID0gYENDUCBtYXBwaW5nYCwgY29sb3VyID0gYENDUCBtYXBwaW5nYCkgKw0KICBnZW9tX2xpbmUoKSArDQogIGdlb21fcG9pbnQoKSArDQogIHNjYWxlX3hfY29udGludW91cyhicmVha3MgPSBzZXEoMTk5MCwyMDIwLDUpKQ0KYGBgDQoNCmBgYHtyfQ0KZ2hnX2VtaXNzaW9uc19kYXRhICU+JSANCiAgZmlsdGVyKEVtaXNzaW9uWWVhciAhPSAiQmFzZVllYXIiKSAlPiUgDQogIG11dGF0ZShFbWlzc2lvblllYXIgPSBhcy5udW1lcmljKEVtaXNzaW9uWWVhcikpICU+JSANCiAgZmlsdGVyKGBOYXRpb25hbCBDb21tdW5pY2F0aW9uIENhdGVnb3JpZXNgICE9IGBDQ1AgbWFwcGluZ2ApICU+JSANCiAgc2VsZWN0KGBOYXRpb25hbCBDb21tdW5pY2F0aW9uIENhdGVnb3JpZXNgLCBgQ0NQIG1hcHBpbmdgKSAlPiUgDQogIHVuaXF1ZSgpDQpgYGANCmNhdGVnb3J5X2lkLGNhdGVnb3J5X25hbWUsc3ViY2F0ZWdvcnlfaWQsc3ViY2F0ZWdvcnlfbmFtZSx5ZWFyLGVtaXNzaW9ucyxlbWlzc2lvbg0KDQpgYGB7cn0NCmVtaXNzaW9uc19zYW5rZXkgPC0gZW1pc3Npb25zX2RhdGEgJT4lIA0KICBzZWxlY3QoY2NwX21hcHBpbmcsIHNvdXJjZV9uYW1lLCBwb2xsdXRhbnQsIGVtaXNzaW9uX3llYXIsIGVtaXNzaW9ucywgdW5pdHMpDQpgYGANCg0KDQpgYGB7cn0NCmdoZ19lbWlzc2lvbnNfY2xlYW4NCmBgYA0KDQoNCmBgYHtyfQ0KZmlsdGVyZWRfZGYgPC0gZ2hnX2VtaXNzaW9uc19jbGVhbiAlPiUgDQogIHNlbGVjdChjY3BfbWFwcGluZywgc291cmNlX25hbWUsIHBvbGx1dGFudCwgZW1pc3Npb25feWVhciwgZW1pc3Npb25zLCB1bml0cykgJT4lIA0KICBmaWx0ZXIocG9sbHV0YW50ID09ICJDTzIiKSAlPiUgDQogIGZpbHRlcihlbWlzc2lvbl95ZWFyID09ICIyMDA1IikNCmBgYA0KDQpgYGB7cn0NCnRvdGFsX2VtaXNzaW9uc19mb3JfZ2FzIDwtIGZpbHRlcmVkX2RmICU+JSANCiAgc3VtbWFyaXNlKHN1bShlbWlzc2lvbnMpKSAlPiUgDQogIHB1bGwoKQ0KDQp0b3RhbF9lbWlzc2lvbnNfYnlfY2F0ZWdvcnkgPC0gZmlsdGVyZWRfZGYgJT4lIA0KICBzZWxlY3QoLXBvbGx1dGFudCkgJT4lIA0KICBncm91cF9ieShjY3BfbWFwcGluZykgJT4lIA0KICBzdW1tYXJpc2UoY2F0X3N1bSA9IHN1bShlbWlzc2lvbnMpLCAuZ3JvdXBzID0gJ2Ryb3BfbGFzdCcpDQpgYGANCg0KYGBge3J9DQpjYXRlZ29yaWVzIDwtIGZpbHRlcmVkX2RmICU+JSANCiAgZGlzdGluY3QoY2NwX21hcHBpbmcpICU+JSANCiAgcHVsbCgpDQoNCm5fY2F0ZWdvcmllcyA8LSBsZW5ndGgoY2F0ZWdvcmllcykNCg0Kc291cmNlcyA8LSBmaWx0ZXJlZF9kZiAlPiUgDQogIGRpc3RpbmN0KHNvdXJjZV9uYW1lKSAlPiUgDQogIHB1bGwoKQ0KDQpuX3NvdXJjZXMgPC0gbGVuZ3RoKHNvdXJjZXMpDQpgYGANCg0KYGBge3J9DQpuX3NvdXJjZXMNCmBgYA0KDQoNCmBgYHtyfQ0Kbm9kZV9uYW1lcyA8LSBjKCJUb3RhbCIsIGNhdGVnb3JpZXMsIHNvdXJjZXMsICJPdGhlciIpDQoNCm5vZGVfbmFtZXNfZGYgPC0gZGF0YS5mcmFtZSgibmFtZSIgPSBub2RlX25hbWVzKQ0KDQp0b3RhbF9zYW5rZXlfdGliYmxlIDwtIHRvdGFsX2VtaXNzaW9uc19ieV9jYXRlZ29yeSAlPiUNCiAgbXV0YXRlKHRvdGFsID0gIlRvdGFsIikgJT4lIA0KICBtdXRhdGUodG90YWwgPSBtYXRjaCh0b3RhbCwgbm9kZV9uYW1lcykgLTEpICU+JSANCiAgbXV0YXRlKGNjcF9tYXBwaW5nID0gbWF0Y2goY2NwX21hcHBpbmcsIG5vZGVfbmFtZXMpIC0xKSAlPiUgDQogIHNlbGVjdChzb3VyY2UgPSB0b3RhbCwNCiAgICAgICAgIHRhcmdldCA9IGNjcF9tYXBwaW5nLA0KICAgICAgICAgdmFsdWUgPSBjYXRfc3VtKQ0KDQp0b3RhbF9maWx0ZXJlZF9lbWlzc2lvbnMgPC0gdG90YWxfc2Fua2V5X3RpYmJsZSAlPiUgDQogIHN1bW1hcmlzZShzdW0odmFsdWUpKSAlPiUgDQogIHB1bGwoKQ0KDQpvdGhlcl9lbWlzc2lvbnMgPC0gdG90YWxfZW1pc3Npb25zX2Zvcl9nYXMgLSB0b3RhbF9maWx0ZXJlZF9lbWlzc2lvbnMNCg0KdG90YWxfb3RoZXJfc2Fua2V5X3RpYmJsZSA8LSB0aWJibGUoDQogICJzb3VyY2UiID0gYygwKSwNCiAgInRhcmdldCIgPSAobWF0Y2goIk90aGVyIiwgbm9kZV9uYW1lcykgLTEpLA0KICAidmFsdWUiID0gYyhvdGhlcl9lbWlzc2lvbnMpDQopDQoNCg0Kc3ViX3NhbmtleV90aWJibGUgPC0gZmlsdGVyZWRfZGYgJT4lIA0KICBzZWxlY3QoLSBjKHVuaXRzLCBwb2xsdXRhbnQsIGVtaXNzaW9uX3llYXIpKSAlPiUgDQogIG11dGF0ZShjY3BfbWFwcGluZyA9IG1hdGNoKGNjcF9tYXBwaW5nLCBub2RlX25hbWVzKSAtMSwNCiAgICAgICAgIHNvdXJjZV9uYW1lID0gbWF0Y2goc291cmNlX25hbWUsIG5vZGVfbmFtZXMpIC0xKQ0KDQpuYW1lcyhzdWJfc2Fua2V5X3RpYmJsZSkgPSBjKCJzb3VyY2UiLCAidGFyZ2V0IiwgInZhbHVlIikNCg0Kc2Fua2V5X3RpYmJsZSA8LSB0b3RhbF9zYW5rZXlfdGliYmxlICU+JSANCiAgYmluZF9yb3dzKHN1Yl9zYW5rZXlfdGliYmxlKSAlPiUgDQogIGJpbmRfcm93cyh0b3RhbF9vdGhlcl9zYW5rZXlfdGliYmxlKQ0KDQpsaW5rc19tYXRyaXggPC0gZGF0YS5mcmFtZShhcy5tYXRyaXgoc2Fua2V5X3RpYmJsZSwgYnlyb3cgPSBUUlVFLCBuY29scyA9IDMpKQ0KDQojIEFkZCBhICdncm91cCcgY29sdW1uIHRvIGVhY2ggY29ubmVjdGlvbjoNCmxpbmtzIDwtIGxpbmtzX21hdHJpeCAlPiUgDQogIG11dGF0ZShncm91cCA9IGNhc2Vfd2hlbigNCiAgICBzb3VyY2UgPT0gMCB+IHBhc3RlKCJ0eXBlXyIsIHRhcmdldCwgc2VwID0gIiIpLA0KICAgIHNvdXJjZSE9MCB+IHBhc3RlKCJ0eXBlXyIsIHNvdXJjZSwgc2VwID0gIiIpDQogICkpDQoNCm5vZGVzIDwtIG5vZGVfbmFtZXNfZGYNCiMgQWRkIGEgJ2dyb3VwJyBjb2x1bW4gdG8gZWFjaCBub2RlLg0KIyBBbGwgb2YgdGhlbSBpbiB0aGUgc2FtZSBncm91cCB0byBtYWtlIHRoZW0gdGhlIHNhbWUgY29sb3VyDQpub2RlcyRncm91cCA8LSBhcy5mYWN0b3IoYygibXlfdW5pcXVlX2dyb3VwIikpDQoNCmVtaXNzaW9ucyA8LSBsaXN0KCkNCg0KZW1pc3Npb25zJG5vZGVzIDwtIG5vZGVzDQplbWlzc2lvbnMkbGlua3MgPC0gbGlua3MNCg0KDQpgYGANCg0KDQoNCg0KDQoNCmBgYHtyfQ0KdG90YWxfZmlsdGVyZWRfZW1pc3Npb25zIDwtIHRvdGFsX3NhbmtleV90aWJibGUgJT4lIA0KICBzdW1tYXJpc2Uoc3VtKHZhbHVlKSkgJT4lIA0KICBwdWxsKCkNCg0Kb3RoZXJfZW1pc3Npb25zIDwtIHRvdGFsX2VtaXNzaW9uc19mb3JfZ2FzIC0gdG90YWxfZmlsdGVyZWRfZW1pc3Npb25zDQoNCnRvdGFsX290aGVyX3NhbmtleV90aWJibGUgPC0gdGliYmxlKA0KICAic291cmNlIiA9IGMoMCksDQogICJ0YXJnZXQiID0gKG1hdGNoKCJPdGhlciIsIG5vZGVfbmFtZXMpIC0xKSwNCiAgInZhbHVlIiA9IGMob3RoZXJfZW1pc3Npb25zKQ0KKQ0KDQpzdWJfc2Fua2V5X3RpYmJsZSA8LSBmaWx0ZXJlZF90aWJibGUgJT4lIA0KICBzZWxlY3QoLWNhdGVnb3J5X2lkLCAtc3ViY2F0ZWdvcnlfaWQsIC15ZWFyKSAlPiUgDQogIG11dGF0ZShjYXRlZ29yeV9uYW1lID0gbWF0Y2goY2F0ZWdvcnlfbmFtZSwgbm9kZV9uYW1lcykgLTEsDQogICAgICAgICBzdWJjYXRlZ29yeV9uYW1lID0gbWF0Y2goc3ViY2F0ZWdvcnlfbmFtZSwgbm9kZV9uYW1lcykgLTEpDQoNCm5hbWVzKHN1Yl9zYW5rZXlfdGliYmxlKSA9IGMoInNvdXJjZSIsICJ0YXJnZXQiLCAidmFsdWUiKQ0KDQpzYW5rZXlfdGliYmxlIDwtIHRvdGFsX3NhbmtleV90aWJibGUgJT4lIA0KICBiaW5kX3Jvd3Moc3ViX3NhbmtleV90aWJibGUpICU+JSANCiAgYmluZF9yb3dzKHRvdGFsX290aGVyX3NhbmtleV90aWJibGUpDQoNCmxpbmtzX21hdHJpeCA8LSBkYXRhLmZyYW1lKGFzLm1hdHJpeChzYW5rZXlfdGliYmxlLCBieXJvdyA9IFRSVUUsIG5jb2xzID0gMykpDQoNCiMgQWRkIGEgJ2dyb3VwJyBjb2x1bW4gdG8gZWFjaCBjb25uZWN0aW9uOg0KbGlua3MgPC0gbGlua3NfbWF0cml4ICU+JSANCiAgbXV0YXRlKGdyb3VwID0gY2FzZV93aGVuKA0KICAgIHNvdXJjZSA9PSAwIH4gcGFzdGUoInR5cGVfIiwgdGFyZ2V0LCBzZXAgPSAiIiksDQogICAgc291cmNlIT0wIH4gcGFzdGUoInR5cGVfIiwgc291cmNlLCBzZXAgPSAiIikNCiAgKSkNCg0Kbm9kZXMgPC0gbm9kZV9uYW1lc19kZg0KIyBBZGQgYSAnZ3JvdXAnIGNvbHVtbiB0byBlYWNoIG5vZGUuDQojIEFsbCBvZiB0aGVtIGluIHRoZSBzYW1lIGdyb3VwIHRvIG1ha2UgdGhlbSB0aGUgc2FtZSBjb2xvdXINCm5vZGVzJGdyb3VwIDwtIGFzLmZhY3RvcihjKCJteV91bmlxdWVfZ3JvdXAiKSkNCg0KZW1pc3Npb25zIDwtIGxpc3QoKQ0KDQplbWlzc2lvbnMkbm9kZXMgPC0gbm9kZXMNCmVtaXNzaW9ucyRsaW5rcyA8LSBsaW5rcw0KYGBgDQoNCg0KDQoNCg0KDQoNCg0KYGBge3J9DQptYWtlX3NhbmtleV9kZnMgPC0gZnVuY3Rpb24oZGF0YSwgdXNlclllYXIsIHVzZXJHYXMpIHsNCiAgbl9jYXRlZ29yaWVzIDwtIGRhdGEgJT4lIA0KICAgIGRpc3RpbmN0KGNhdGVnb3J5X25hbWUpICU+JSANCiAgICBucm93KCkNCiAgDQogIHRvdGFsX2VtaXNzaW9uc19mb3JfZ2FzIDwtIGRhdGEgJT4lIA0KICAgIGZpbHRlcihlbWlzc2lvbiA9PSB1c2VyR2FzKCkpICU+JSANCiAgICBmaWx0ZXIoeWVhciA9PSB1c2VyWWVhcigpKSAlPiUNCiAgICBzdW1tYXJpc2Uoc3VtKGVtaXNzaW9ucykpICU+JSANCiAgICBwdWxsKCkNCiAgDQogIGZpbHRlcmVkX3RpYmJsZSA8LSBkYXRhICU+JQ0KICAgIGZpbHRlcihlbWlzc2lvbiA9PSB1c2VyR2FzKCkpICU+JSANCiAgICBzZWxlY3QoLWVtaXNzaW9uKSAlPiUgDQogICAgZmlsdGVyKHllYXIgPT0gdXNlclllYXIoKSkgJT4lIA0KICAgIGZpbHRlcihlbWlzc2lvbnMgPiB1c2VyUmVzb2x1dGlvbigpKQ0KICANCiAgdG90YWxfZW1pc3Npb25zX2J5X2NhdCA8LSBmaWx0ZXJlZF90aWJibGUgJT4lDQogICAgZ3JvdXBfYnkoY2F0ZWdvcnlfbmFtZSkgJT4lIA0KICAgIHN1bW1hcmlzZShjYXRfc3VtID0gc3VtKGVtaXNzaW9ucyksIC5ncm91cHMgPSAnZHJvcF9sYXN0JykNCiAgDQogIGNhdGVnb3JpZXMgPC0gZmlsdGVyZWRfdGliYmxlICU+JQ0KICAgIGRpc3RpbmN0KGNhdGVnb3J5X25hbWUpICU+JSANCiAgICBwdWxsKCkNCiAgDQogIHN1YmNhdGVnb3JpZXMgPC0gZmlsdGVyZWRfdGliYmxlICU+JQ0KICAgIGRpc3RpbmN0KHN1YmNhdGVnb3J5X25hbWUpICU+JSANCiAgICBwdWxsKCkNCiAgDQogIG5vZGVfbmFtZXMgPC0gYygiVG90YWwiLCBjYXRlZ29yaWVzLCBzdWJjYXRlZ29yaWVzLCAiT3RoZXIiKQ0KICANCiAgbm9kZV9uYW1lc19kZiA8LSBkYXRhLmZyYW1lKCJuYW1lIiA9IG5vZGVfbmFtZXMpDQogIA0KICB0b3RhbF9zYW5rZXlfdGliYmxlIDwtIHRvdGFsX2VtaXNzaW9uc19ieV9jYXQgJT4lDQogICAgbXV0YXRlKHRvdGFsID0gIlRvdGFsIikgJT4lIA0KICAgIG11dGF0ZSh0b3RhbCA9IG1hdGNoKHRvdGFsLCBub2RlX25hbWVzKSAtMSkgJT4lIA0KICAgIG11dGF0ZShjYXRlZ29yeV9uYW1lID0gbWF0Y2goY2F0ZWdvcnlfbmFtZSwgbm9kZV9uYW1lcykgLTEpICU+JSANCiAgICBzZWxlY3Qoc291cmNlID0gdG90YWwsDQogICAgICAgICAgIHRhcmdldCA9IGNhdGVnb3J5X25hbWUsDQogICAgICAgICAgIHZhbHVlID0gY2F0X3N1bSkNCiAgDQogIHRvdGFsX2ZpbHRlcmVkX2VtaXNzaW9ucyA8LSB0b3RhbF9zYW5rZXlfdGliYmxlICU+JSANCiAgICBzdW1tYXJpc2Uoc3VtKHZhbHVlKSkgJT4lIA0KICAgIHB1bGwoKQ0KICANCiAgb3RoZXJfZW1pc3Npb25zIDwtIHRvdGFsX2VtaXNzaW9uc19mb3JfZ2FzIC0gdG90YWxfZmlsdGVyZWRfZW1pc3Npb25zDQogIA0KICB0b3RhbF9vdGhlcl9zYW5rZXlfdGliYmxlIDwtIHRpYmJsZSgNCiAgICAic291cmNlIiA9IGMoMCksDQogICAgInRhcmdldCIgPSAobWF0Y2goIk90aGVyIiwgbm9kZV9uYW1lcykgLTEpLA0KICAgICJ2YWx1ZSIgPSBjKG90aGVyX2VtaXNzaW9ucykNCiAgKQ0KICANCiAgc3ViX3NhbmtleV90aWJibGUgPC0gZmlsdGVyZWRfdGliYmxlICU+JSANCiAgICBzZWxlY3QoLWNhdGVnb3J5X2lkLCAtc3ViY2F0ZWdvcnlfaWQsIC15ZWFyKSAlPiUgDQogICAgbXV0YXRlKGNhdGVnb3J5X25hbWUgPSBtYXRjaChjYXRlZ29yeV9uYW1lLCBub2RlX25hbWVzKSAtMSwNCiAgICAgICAgICAgc3ViY2F0ZWdvcnlfbmFtZSA9IG1hdGNoKHN1YmNhdGVnb3J5X25hbWUsIG5vZGVfbmFtZXMpIC0xKQ0KICANCiAgbmFtZXMoc3ViX3NhbmtleV90aWJibGUpID0gYygic291cmNlIiwgInRhcmdldCIsICJ2YWx1ZSIpDQogIA0KICBzYW5rZXlfdGliYmxlIDwtIHRvdGFsX3NhbmtleV90aWJibGUgJT4lIA0KICAgIGJpbmRfcm93cyhzdWJfc2Fua2V5X3RpYmJsZSkgJT4lIA0KICAgIGJpbmRfcm93cyh0b3RhbF9vdGhlcl9zYW5rZXlfdGliYmxlKQ0KICANCiAgbGlua3NfbWF0cml4IDwtIGRhdGEuZnJhbWUoYXMubWF0cml4KHNhbmtleV90aWJibGUsIGJ5cm93ID0gVFJVRSwgbmNvbHMgPSAzKSkNCiAgDQogICMgQWRkIGEgJ2dyb3VwJyBjb2x1bW4gdG8gZWFjaCBjb25uZWN0aW9uOg0KICBsaW5rcyA8LSBsaW5rc19tYXRyaXggJT4lIA0KICAgIG11dGF0ZShncm91cCA9IGNhc2Vfd2hlbigNCiAgICAgIHNvdXJjZSA9PSAwIH4gcGFzdGUoInR5cGVfIiwgdGFyZ2V0LCBzZXAgPSAiIiksDQogICAgICBzb3VyY2UhPTAgfiBwYXN0ZSgidHlwZV8iLCBzb3VyY2UsIHNlcCA9ICIiKQ0KICAgICkpDQogIA0KICBub2RlcyA8LSBub2RlX25hbWVzX2RmDQogICMgQWRkIGEgJ2dyb3VwJyBjb2x1bW4gdG8gZWFjaCBub2RlLg0KICAjIEFsbCBvZiB0aGVtIGluIHRoZSBzYW1lIGdyb3VwIHRvIG1ha2UgdGhlbSB0aGUgc2FtZSBjb2xvdXINCiAgbm9kZXMkZ3JvdXAgPC0gYXMuZmFjdG9yKGMoIm15X3VuaXF1ZV9ncm91cCIpKQ0KICANCiAgZW1pc3Npb25zIDwtIGxpc3QoKQ0KICANCiAgZW1pc3Npb25zJG5vZGVzIDwtIG5vZGVzDQogIGVtaXNzaW9ucyRsaW5rcyA8LSBsaW5rcw0KICANCiAgcmV0dXJuKGVtaXNzaW9ucykNCn0NCmBgYA0KDQoNCmBgYHtyfQ0KbWFrZV9zYW5rZXlfZGZzKGVtaXNzaW9uc19zYW5rZXksIHVzZXJZZWFyID0gMjAwNSwgdXNlckdhcyA9ICJDSDQiLCB1c2VyUmVzb2x1dGlvbiA9IDUwKQ0KYGBgDQpgYGB7cn0NCg0KYGBgDQoNCg==